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Abstract 

§ . Recent lattice data from CP-PACS, UKQCD, SESAM/TXL and the Pisa group 

CN , regarding the quark mass dependence of the topological susceptibility in 2-flavour 

I QCD are compared to each other and to theoretical expectations. The latter get 

specified by referring to the QCD finite-volume partition function with "granularity" 
. which accounts for the entropy brought by instantons and anti-instantons. The chiral 

! condensate in Nf = 2 QCD, if determined by this method, turns out surprisingly large. 
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_ 1 Introduction 
o 

O '. Quenched QCD is, by many respects, a reasonable approximation to real ("full") QCD; in 
particular for the light meson spectrum this approximation was found to work so well that 
a major effort had to be made to drive statistical errors small enough so that significant 
■ differences between the computed quenched and the experimentally observed (unquenched) 
^ . spectrum could be seen On the other hand, there are observables which prove highly 
D ! sensitive on quenching effects. An important example of this latter category is the topolo- 
^ ' gical susceptibility which, in the continuum, is defined^, through [§] 

>< ■ X= fd'z dMT{k^{x)k%y)}) , (1) 

where z = x—y and 

k'ix) = J^e'^-M:(a;) {d,Al{x) - |r^^Aj(x) (2) 
denotes the Chern current which is related to the topological charge density via 



g(j^) = a^fli) = J^G;„(i)G""'(i) with G"'' = ie'"'"''G,, . (3) 



The topological susceptibility x plays a prominent role in QCD, because its nonvanishing 
quenched counterpart Xoo explains through the Witten-Veneziano relation 

2j\j 

Ml, - (2Mi - Ml) = -^xoo (4) 

^Upon including a factor e"^^ into the definition x would pick up an additional dependence on . 
hi this article we shall only be interested in x ^t zero virtuality. 
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to leading order in 1/Nc why the t]' is heavier than the octet of pseudoscalar mesons and 
hence not a pseudo-Goldstone. The topological susceptibility x iii the full theory is an 
interesting object to study unquenching effects, because it may deviate from the quenched 
prediction Xoo by up to 100%, if the dynamical quarks tend to be light. In this respect 
one should note that the dependence of x = x{^fy''^) the number of dynamical quarks 
and their masses is not explicit in the definition (||); it enters implicitly via the vacuum 
expectation value (...) which includes the weight of the fermion functional determinant. 
This means that the topological susceptibility is useful to study the QCD vacuum structure 
and how the latter is affected by the presence of light dynamical quarks. It is therefore little 
surprise that the first reliable lattice data for the topological susceptibility in full QCD were 
eagerly awaited for. 

Below we shall analyze the data by CP-PACS |], UKQCD 0, SESAM/TXL § and 
the Pisa group @] . We will see that they are roughly consistent with each other, but follow 
theoretical expectations to a limited extent only. Potential sources of error are analyzed, 
and it is shown that specifically for the topological susceptibility finite- volume effects may 
be quantitatively assessed by deriving it from the QCD finite-volume partition function for 
which an elaborate version is given. On the practical side, a non-standard determination of 
the chiral condensate in Nf = 2 QCD is presented, and the result is found surprisingly large. 



2 Theoretical expectations 

To begin, we shall remind ourselves of the theoretical expectations regarding the dependence 
of the topological susceptibility x on the sea-quark masses mi, . . . ,m]\f^.. For simplicity the 
latter shall be taken degenerate, i.e. we are interested in the function x = xi^fi''^) with 
Nf > 2. In the present section the four-volume V of the box is assumed "infinite"; finite- 
volume effects will be analyzed below. 

Detailed predictions on x = xi^fi ^) i^^^J be obtained for the two regimes where the 
quark masses tend to be either very light or very heavy. In the former case the flavour 
singlet axial Ward-Takahashi (WT) identity dictates that x raises linearli^ 0; 0; 0; viz. 

X = ^m + o{m) , (5) 

where S = — limm^o liniy_»oo (V'^) is the chiral (one- flavour) condensate in the chiral limit[| 
and the o{m) term contains higher order corrections^ which may be worked out explicitly 
in chiral perturbation theory |]T0|. On the other hand, for heavy dynamical quarks the 



topological susceptibility curve will be almost flat, because it must gradually approach the 
quenched value Xoo which is finite. A numerical estimate of the latter quantity by means of 
the Witten-Veneziano relation (|) yields Xoo ^ (180 MeV)^ ~ 1.05 lO^MeV^. In addition, it 
is easy to come up with plausible arguments that x grows monotonically as the quark mass 
increases and that the second derivative is negative'^. 

Hence, without doing a detailed analysis, a clear picture emerges for the qualitative 
behaviour of the topological susceptibility as a function of the quark mass over the whole 
range of the latter. Even the transition regime may be estimated by just equating the 

^Here the assumption V ^oo proves essential; for the behaviour at finite V see below. 

•^Both E and m are scheme- and scale-dependent; the product, however, is an RG-invariant quantity. In 
this article, individual quotes for E, m mean that standard conventions (MS, /x = 2 GeV) have been adopted. 

■^The reason why I write o{m) rather than Olm?) is that the higher order terms may contain chiral 
logarithms; a chiral expansion is (in general) non-analytic. 
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Figure 1: Basic illustrations regarding the topological susceptibility [in MeV^j versus the quark 
mass [in MeV] for 2-fiavour QCD with S = (288MeV)3 (cf. footnote 5). LHS: The leading order 
chiral prediction in the large Leutwyler-Smilga regime without instanton effects (jsj) (raising slope) 
and the quenched prediction from the Witten-Veneziano relation (Q) (Eat line). RHS: The curves 
with n = 1,2,4, 8 (from the lowest to the highest) out of the phenomenological family (^), all of 
which share the required asymptotic behaviour. 



leading-order chiral expression @ to the asymptotic value Xoo, as is illustrated on the l.h.s. 
of Fig. 1: Upon usingQ S = (288 MeV)^ = 2.4 10^ MeV^ one gets mtrans ^ 88 MeV for 2-flavour 
QCD, which is accessible in current dynamical simulations. 

What these considerations cannot yield is the functional form of the topological sus- 
ceptibility curve. Of course it is easy to come up with a whole bunch of phenomenological 
models or ansatze which interpolate in various ways between the asymptotic regimes; e.g. 

x{Nf, m) = ((S m/Nf)- + Xoo")"'/" (6) 

is a one-parameter family of interpolating curves (some of which are displayed in the r.h.s. 
of Fig. 1) with the required asymptotic behaviour. 

In summary, one would expect up-to-date simulations of full 2-fiavour QCD to show the 
turn-over from a linear to an almost fiat behaviour at a quarlc mass of 0(88 MeV), and to 
select the most promising phenomenological overall fitting curve from the set (j^). 



3 Lattice data 

Before considering the actual data, we should remind ourselves that it is not the original 
definition (|ID of the topological susceptibility which is typically used on the lattice, but a 
simplified version which - at least in the continuum - proves equivalent. 

The simplification talces place as the result of a two-step procedure. The first step is to 
assume translation invariance of the vacuum, so (^ simplifies to 

X = Jd'x{T{q{x)q{0)}) , (7) 

where we have also assumed that the (5-functions due to the derivatives in (|1|) hitting the 
T-ordering do not contribute after integration. The second step assumes that there is, in 
addition, no correlation between topological charge densities, so (0) simplifies to 

X^^, (8) 

^The numerical value stems from the Gell-Mann-Oakes-Renner (GOR) PCAC-relation M^F^ ~ 2Y,m, 
together with 139 MeV, ~93 MeV and m~ 3.5 MeV (cf. footnote 3) [|lll. 
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Figure 2: CP-PACS data for the distribution of topological charges on a 16^ x 32 lattice after 
cooling with an RG-improved action. The configurations were generated using the same (Iwasaki) 
action at fixed /3 = 1.95 together with 2 Eavours of mean-field improved clover quarks (c^^ = 1.53) 
at 4 different n-values. Figure by courtesy of the CP-PACS collaboration. 



where 



327r2 



d^xG^'^^G'"'^ (9) 



is the (global) topological charge which, in the continuum, takes an integer value. The 
important point is that the original definition (||) provides a detailed prescription of how 
to deal with the (5-like singularity in the product q{x)q{0) which is not explicit any more in 
the simplified version (Q). As a matter of consequence []12[, a lattice implementation of (|1]) 
is subject to a multiplicative renormalization only 

X = Z'a^X + Oia'') (10) 

with the renormalization factor of the topological charge 

q = Zq + Oia'), (11) 

whereas a lattice version of (|^) shows an additional additive renormalization, viz. 

X = M{P, Nf, m) + Z^a^x + 0(0^) . (12) 

While in practice both the multiplicative renormalization factor Z and the additive piece M 
do not cause any difficulties, since they may be determined by means of heating techniques 
T3| , an explicit verification of the three definitions (|l], 0, P) proving equivalent after they 



have been implemented on the lattice is still being awaited. Since, as we have seen, proper 
treatment of contact terms is essential in this game (cf. |]14| and references therein), such 
an investigation might be interesting in its own right. 

After these cautionary remarks we are ready to esteem the latest lattice data which 
use the simplified definition (|^), thereby dealing - in principle - with the renormalization 
pattern (|T^). In order to get a feeling of what is actually done, it is useful to consider 
the distribution of topological charges as found by the CP-PACS collaboration 0|. Fig. 
2 shows the histograms they have gotten on a 16^ x 32 lattice at /? = 1.95 (using the 
Iwasaki action) and at four different K-values (/t = 0.1375,0.1390,0.1400,0.1410, giving 
M^/Mp = 0.80,0.75,0.69,0.58, respectively) for 2 flavours of mean-fleld improved clover 
quarks (Cg^ = 1.53) after cooling with the Iwasaki action and binning q to integer values. 
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Figure 3: CP-PACS data for the topological susceptibility [in units of r^^] versus the quark mass 
[expressed in terms of M^Tq] in 2-fiavour QCD, as determined from the q-distributions shown in 
Fig. 2. The linearly raising slope indicates the asymptotic behaviour for light dynamical quarks 
(Xoro ~ {F'^rl/2Nf)Mlrl ~ 0.0133 M^rg for ~ 93 MeV and tq ~ 0.49 fm ), and the horizontal line 
the one for (infinitely) heavy dynamical quarks (Xoo'''o — (180 MeV)^ (0.49 fm)'' ~ 0.04j, as in the 
l.h.s. of Fig. 1. The data are expected to lie in the sector beneath either one of these lines, i.e. 
XrQ<min(xo^O'/too^o)' w-iiere Xoo might be somewhat larger (cf. text). Data taken from 



Since the cooling procedure brings both renormahzation constants [Z and M) close to 
their continuum values (1 and 0, respectively), one can basically read off the expectation 
values (g^) from the widths of the distributions shown in Fig. 2. This holds, because 

for a gaussian distribution. Hence, from a glimpse at Fig. 2 one might get the impression 
that the topological susceptibility gets reduced as the quark mass decreases, since the q- 
distribution does indeed get narrower as k increases. What this naive consideration neglects 
is that the topological susceptibility x is defined as the variance of these gaussian curves 
divided by the physical four- volume of the box (cf. (||, |13D), and the latter changes with 



K, because in full QCD the physical lattice spacing a is a non-trivial function of f3 and k. 
To be specific: a decreases if k gets increased at fixed /3, and for the volume the reduction 
is particularly strong, as it involves the fourth power of a. Taking this effect into account 
yields the data for x versus the dynamical quark mass m in = 2 QCD shown in Fig. 3. 
They seem to indicate an almost flat topological susceptibility curve, and this means that 
the narrowing of the g-distribution in Fig. 2 was exclusively due to the volume V decreasing 
with increasing k (at fixed (3) and not due to a decrease of x- Iii other words: The data in 
Fig. 3 show no sign of the expected transition to the linear behaviour (^) characteristic for 
light dynamical quarks, even though the most chiral point is at M^rl = 1.6Q, i.e. at a quark 
mass of m~50MeV (using = 93 MeV, S = (288 MeV)^) which is well below our estimate 
of the transition regime; we had mtrans — 88MeV for Nf = 2. The troublesome point is the 
most chiral one; it seems to violate the constraint x^Xo with xo defined as the r.h.s. of (^. 
Potential reasons for the supposed excess of the measured x will be discussed below. 

From Fig. 4 one sees that the data generated by the CP-PACS collaboration on the 
larger lattice (24^x48 at (3 = 2.10) as well as the results obtained by UKQCD and the Pisa 
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Figure 5: Collection of recently published data for the topological susceptibility in Nf = 2 QCD. 



group are "better" in the sense that they barely support an entirely flat behaviour, but they 
still seem to be at the edge of being compatible with the constraints imposed by the known 
asymptotic behaviour for small quark masses. For the CP-PACS data at (3 = 2.10 and the 
UKQCD data the hypothesis that the true 2-fiavour topological susceptibility curve lies in 
the sector beneath either one of the two indicated lines seems conceivable, in particular if 
one takes into account that the horizontal line itself is subject to some uncertainty^. The 
remaining constraints - the first derivative of x(m) must always be positive, the second 
always negative - seem to be obeyed with the (unsignificant) exception of one point in 
each one of these two simulations. Thanks to relatively large error bars, it seems that 
the data obtained by the SESAM-coUaboration are at the edge of being compatible with 
the two asymptotic constraints. For the data due to the Pisa group the situation looks 
less promising: Each one of their data points is compatible, on a 2o"-level, with the linear 
constraint applicable in the deep chiral regime, but the combination of the five points is 
not. For obvious reasons nothing can be said about the first two derivatives from either the 
SESAM or the Pisa data, the error bars are just too large. 

Another simple cross-check is to combine the results obtained by the various collabora- 
tions in a single plot and to see whether there is an obvious trend for the data to depend 
on the group, an endeavour undertaken in Fig. 5. The first impression is that the data are 
vaguely consistent with each other, but even for the combined data it is not clear whether 
they show all of the transient behaviour we expect to see with active flavours which are 
neither light enough to produce a linear curve through the origin nor heavy enough to have 
the topological susceptibility flatten out at the saturated quenched level: They seem to 
indicate a positive slope (though the band is rather fuzzy), but it is not clear whether the 
combined data do exhibit the negative curvature genuine to the transition region anticipated 
at m = 0(88MeV) or M^r^ = 0{3) for Nf = 2 (cf. Fig. 1). 

The consistency between the various data sets cannot be discussed, however, without 
taking into account that the groups have adopted rather different simulation strategies. 



^The dotted line assumes Xoo = (190 MeV)''. Using Xoo = (180 MeV)^ or Xoo = (200 MeV)"* shifts the Une 
to 0.04 or 0.06, respectively. Note that the situation for the chiral constraint is different; here additional 
positive terms on the r.h.s. of the GOR-relation M^F^ = 2Em + ... would reduce its slope, and the bound 
on S is relatively tight, since the experimental determination F^^c^idSMeY is rather precise. 



7 



CP-PACS-1 


16^x32 /3=1.95 Cg^=1.53 


K 


0.1410 0.1400 0.1390 0.1375 
0.586 0.688 0.751 0.805 


Mlrl 

TT U 

4 

a„[im\ 
x = VT,m 


1.66 2.83 3.74 5.00 
0.041 0.038 0.039 0.044 
0.170 0.181 0.193 0.204 

109. 141. 182. 227. 

42.0 92.1 157. 263. 



CP-PACS-2 


24^x48 p = 2.10 Cs^=1.47 


MJM, 


0.1382 0.1374 0.1367 0.1357 
0.575 0.690 0.757 0.806 


Mlrl 

4 

X^Q 

a^[fm] 
x = VIlm 


1.76 3.23 4.44 5.86 
0.032 0.047 0.057 0.045 
0.113 0.120 0.126 0.134 

107. 138. 167. 215. 

43.7 103. 171. 292. 



UKQCD 


W X 32 


partially 


matched 








5.20 


5.20 


5.20 


5.26 


5.29 




2.02 


2.02 


2.02 


1.95 


1.92 


K 


0.13565 


0.1355 


0.1350 


0.1345 


0.1340 


MjMp 


0.567 


0.600 


0.688 


0.791 


0.830 


MM 

TV U 


1.90 


1.98 


3.49 


5.49 


6.61 


xrt 


0.023 


0.031 


0.040 


0.043 


0.064 




0.094 


0.097 


0.103 


0.104 


0.102 




10.3 


11.7 


14.8 


15.4 


14.1 


x = V'Lm 


9.1 


10.7 


23.9 


39.1 


43.1 



SESAM 


16^x32 p = 5.6 Csw=0 


243 X 40 dito 


MJM, 


0.1575 0.1570 0.1565 0.1560 
0.692 0.763 0.813 0.834 


0.1580 0.1575 
0.574 0.704 


Mlrl 

xrt 
a-ro [fm] 
V[fm^] 
x = VIlm 


2.81 3.56 4.45 5.20 
0.050 0.051 0.065 0.047 
0.084 0.091 0.095 0.098 
6.5 9.1 10.6 12.1 
8.5 15.1 21.8 29.1 


1.53 2.66 
0.047 0.053 
0.080 0.085 
23.1 28.9 
16.4 35.6 



PISA 


16^ am = 0.0125 (iV/ = 2 stag.) 


(3 


5.40 5.50 5.55 5.60 5.70 


Mlrl 

xrt 
aa [fm] 

x = VIlm 


0.54 0.69 0.77 0.89 1.12 
0.024 0.032 0.023 0.038 0.028 
0.174 0.134 0.125 0.108 0.086 
60.3 21.3 15.9 8.78 3.55 
15.1 6.8 5.7 3.6 1.8 



Tables 1-5: Key data of the five studies. For error bars see Figs. 1-5 and Refs. jj^, |^ ^, In all 
cases the Leutwyler-Smilga parameter x = VSm has been estimated by re-expressing it, through 
the GOR relation, as x = V MIF"^ /2 = V{MItI){F^ /2rl) with F^~93MeV,ro~0.49fm. 
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CP-PACS has decided to simulate at fixed (3, going more chiral by tuning k closer to Kcrit- 
As discussed above, this approach reduces the physical lattice spacing a for low M^^jMp. In 
order to keep a constant, one would have to reduce /3 appropriately while increasing k, an 
approach which - if successfully implemented - produces "matched ensembles" . This is the 
philosophy adopted by UKQCD. For technical reasons (Cg^ is not known nonperturbatively 
for /? < 5.2) this strategy had to be sacrificed for their two most chiral points, i.e. their 
data set is, strictly speaking, just "partially matched". SESAM, on the other hand, has 
decided to take everywhere the simplest choice, i.e. to vary k at fixed /3 while not making 
any attempt to improve the Wilson fermions. An orthogonal strategy has been adopted by 
the Pisa group: They simulate (using the staggered rather than a Wilson-type formulation 
for their dynamical flavours) at fixed am, thus decreasing (5 in order to go more chiral. 
Note that - contrary to what happens in a series of fixed-/? simulations - the Pisa approach 
makes the lattice coarser as one goes more chiral. 

A possible inconsistency between the results by the different groups could come from a 
variety of reasons: 

(a) Scaling violations: Because of the computational costs involved in full QCD simu- 
lations the latter are typically performed on relatively coarse lattices (cf. Tables 1-5). For 
the simulations analyzed in this article scaling violations are supposed to be "small" since 
0(a) effects have been removed by mean-field techniques (CP-PACS) or non-perturbatively 
(UKQCD) or they might be numerically suppressed due to a small a (SESAM) or absent 
anyway (Pisa group). Nontheless, the physical lattice spacing varying from one data point 
to another in the same study might imply unequal shifts due to cutoff effects. 

(b) Chirality violation with Wilson-type fermions: Wilson fermions (SESAM) violate 
the chiral symmetry at nonzero lattice spacing, and this holds true, too, for all ultralocally 
improved descendents (e.g. the clover-fermions used by CP-PACS and UKQCD). 

(c) Finite-volume effects: From Tables 1-5 one learns that the five studies have been 
conducted in boxes of rather different physical four-volume V and at rather different values 
of the Leutwyler-Smilga parameter x = VT,m. Hence, if finite- volume effects are not fully 
under control, potential inconsistencies of type (a) or (b) might get masked by them. 

I have a few comments on each one of these points: 

On (a): UKQCD claimed that in order to get simulation results which indicate a quark 
mass dependence of the topological susceptibility it was crucial to use the strategy of 
generating matched ensembles, rather than running at fixed (3 0] . Despite this assertion it 
holds true that in principle either method should work: If one could do both a continuum 
and a chiral extrapolation, it should not matter at all where in the (/3, k) plane the data were 
collected. In practice a continuum extrapolation at fixed quark mass is unaffordable at the 
present time, and hence each group has to live with the single data point for a given value 
of M2 rQ. Since both UKQCD and CP-PACS work with 0(a) improved actions for which a 
continuum extrapolation was shown to bring just a minor shift in the case of pseudoscalar 
masses , it is natural to expect the topological susceptibility data to be reasonably close 
to their continuum values too. If unequal magnitudes of scaling violations in a fixed-/? 
strategy were indeed responsible for the majority of the CP-PACS-1 "flatness problem", 
then one would expect the data generated by CP-PACS at /?= 1.95 to be subject to larger 
shifts for large quark masses (where the lattice is coarser) versus more trustworthy at small 
quark masses (where a is smaller), but from Fig. 3 one gets the impression that the problem 
with that particular set is primarily in its most chiral point. Scaling violations may indeed 
account for the bulk of the difference between the UKQCD and CP-PACS-1 data sets, but 
it seems that such a difference is much more likely associated with the respective overall 
differences between typical lattice spacings (Tables 1 and 3 indicate that UKQCD works 
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on much finer lattices than those employed by CP-PACS a.t (3= 1.95) rather than whether 
these sets are matched or not. This is also the conclusion one is led to when comparing 
the CP-PACS-1 data (Table 1) to those entitled CP-PACS-2 (Table 2): The latter seem to 
indicate a quark mass dependence even though they were generated at fixed P = 2.10. Since 
the volume parameters V and x are in one-to-one analogy between CP-PACS-1 and CP- 
PACS-2 (see Tables 1 and 2), the only difference is the lattice spacing which is significantly 
smaller at the higher /3-value. 

On (b): The quark mass dependence of the topological susceptibility in the extreme 
chiral regime is derived from the flavour singlet axial WT-identityQ. With Wilson quarks 
(or descendents which do not satisfy the Ginsparg- Wilson relation) this identity receives 
additional terms [^. This may result in the lattice topological susceptibility x being 
enhanced compared to its continuum counterpart x- It should be mentioned, however, that 
this effect is, in principle, accounted for by the additive renormalization term in (|T2|). Hence, 
on a practical level, an accurate determination of M = M{(3, Nf,m) should be sufficient to 
deal with this challenge. 

On (c): As indicated in footnote 2, for the topological susceptibility to raise linearly 
as a function of the (small) quark mass it is important that the volume is "inflnite", i.e. 
in more technical terms that the Leutwyler-Smilga parameter (|T^) is large. How "large" 
it has to be for this behaviour to show up and what pattern the topological susceptibility 
follows if this criterion is not satisfled will be discussed in the following two sections. As we 
shall see, flnite- volume effects are the only source of contamination of the lattice topological 
susceptibility for which an a-priory quantitative assessment may be achieved. 

4 Leutwyler-Smilga regimes 

Considering a theory which exhibits the spontaneous breakdown of a continuous global 
symmetry (e.g. of the SU{Nf)A group for QCD) in a flnite volume, one expects to see both 
symmetry restoration phenomena and the onset of Goldstone boson production, if the box 
is taken sufficiently "small" or "large", respectively. An obvious question is what sets the 
scale, i.e. by which standards does the box have to be small or large to trigger one or the 
other type of phenomena ? 

The naive guess is the lightest particle at hand, i.e. the mass of the pseudo- Goldstone 
boson produced in the infinite- volume limit. From this one would expect symmetry restora- 
tion to take place for M^L<^1 and SSB to become manifest for Mj^L^l. 

As Leutwyler and Smilga (LS) have shown |p, there is, in addition, the parameter 

X = VSm (14) 

which also depends on the quark mass and which proves particularly useful in a lattice 
context, since it decides which one of the limits m ^ 0,V ^ oo is in the "inner" posi- 
tion and therefore winning^. In other words: For x <^ 1 the chiral symmetry is almost 
restored, i.e. quarks and gluons are the dominant degrees of freedom, whereas for x ^ 1 
the chiral symmetry is (though the box volume is still finite) effectively broken, i.e. long 

''See e.g. Ref. It may also be derived from chiral perturbation theory, since the latter is just a neat 
way of implementing the constraints imposed by the axial WT-identity. 

^In QCD the two limits limj„^Oj hmy^oo are known not to commute: The chiral condensate vanishes if 
m^O is performed first, but it tends to — E, if the volume is sent to infinity first and then the chiral limit 
is taken. In a numerical analysis the two limits cannot be considered separately; the upshot is that lattice 
data will reflect the former situation for x<Cl and the latter for x^l. 
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range Green's functions are dominated by Goldstone excitations. Furthermore, the value of 
X has a bearing on whether standard^ physical observables depend on the net topological 
charge q of the gauge background or not |^: In the small LS-regime (x^ 1) all observables 
depend massively on |g|, and the finite- volume partition function is entirely dominated by 
the contribution from the charge zero sector viz. 

-(5,0 (x<l). (15) 

On the other hand, in the large LS-regime {x ^ 1) standard observables prove basically 



independent of |g| (for a numerical check and some qualifications see [0]), and the finite- 
volume partition function gets broad and gaussian viz. 

It is important to note that the LS-parameter x and the "naive" parameter 

y = M^L (17) 

are not the same: In either one of the regimes with pronounced (x^l) or mild (x~l) SSB 
the box may be large {y ^ 1), intermediate {y ~ 1) or even small {y <^ 1) w.r.t. the pion 
correlation length classification, whereas in the symmetry restoration regime (x^l) such 
a distinction does not make sense, since there the pion is not a useful degree of freedom. 
That X and y cover different aspects of the box being "small" or "large" is also borne out by 
the fact that the minimum box length scales differently as a function of the inverse quark 
mass, if one wants the box to be large w.r.t. either criterion, i.e. 

Lrain OC (1/m)^/^ for m ^ at fixed (large) x (18) 
Lmin oc (1/m)^/^ for m ^ at fixed (large) y, (19) 

which simply follows from the definition ( [I^ ) of x and the GOR relation, respectively. 
Furthermore, it is worth noticing that in a lattice context the LS-classification refers to the 
mass of the sea-quarks (i.e. the dynamical quarks which influence the weight of a gauge 
configuration through the determinant), whereas the conventional one refers to the mass of 
the current-quarks (i.e. those from which hadronic observables are built). 

In QCD the chiral symmetry restoration transition is intricately related to the decon- 
finement transition, and this is the reason why even the combination of the parameters x 
and y is not sufficient to describe the broad properties of the finite-volume system. This 
is easily seen from the fact that the confinement- deconfinement transition takes place at a 
critical temperature Tc~ 275 MeV in pure S'[/(3)-gluodynamics and at a somewhat reduced 
temperature Tc~ 175 MeV in QCD with Nf = 2 massless flavours. In other words: Tc stays 
at the same order of magnitude while both x (cf. (|T^) and y (cf. |T^) decrease from infinity 
to zero as the quarks are tuned from infinitely heavy to weightless. Hence a third scale must 
be involved which - contrary to x and y - barely depends on Nf and m. It is tempting^ to 



^The meaning is: an observable not related to the [/(l)yi-issue, e.g. Af^r, F.^, Mp, Vqq, but not M^'. 
^''This is not intended to overrule the more traditional view that the transition occurs when the hadrons 
are packed so densely that they start to overlap, though it does have the decisive advantage over the latter 
that it does not try to predict the transition from the properties of one phase only. (The true condition 
is of course the equalness of the free energy densities.) In addition, it should be stated that no claim is 
made that instantons are truly responsible for the confinement-deconfinement transition (at least random 
superpositions of instantons do not reproduce the string tension ) , but that - because of the relationship 
the chiral and the confinement transition have in QCD - a complete assessment of whether the box is "small" 
or "large" needs a scale which does not depend on the quark masses. The scale Lmst fulfills this requirement 
and turns out to be numerically in the right order, and that is all we need for the present analysis. 
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identify the missing scale with the one due to instantons, i.e. to introduce the parameter 



z = ^ (20) 

-^inst 

where Lmin is the smallest[3 hnear dimension of the box, and Linst is a typical linear size 
associated to an instanton, e.g. Linst = 2pinst — 0.7fm or Linst = V^nst — ^ [05 since such a 
scale is, to a large extent, independent of the number and mass of the active flavours. For 
z-Cl the system is in the deconfined state, i.e. quarks and gluons represent the appropriate 
degrees of freedom, whereas for z ^ 1 confinement is manifest, i.e. asymptotic states are 
colour singlets, and the associate composite fields are the building blocks of the effective 
description appropriate in that regime. 

As we shall see in the following section, the subtle interplay between the LS-scale x 
(which depends on the quark masses) and the instanton size (which barely does) is crucial 
for an understanding of the overall behaviour of the topological susceptibility curve in QCD. 



5 Finite-volume partition function with "granularity" 

We are now in a position to exploit the fact that one may actually derive the topological 
susceptibility from the finite-volume partition function, and for the latter several versions 
are available in the literature. 



In the large LS-regime the situation is particularly simple, since from the representation 
(p!6|) of the QCD finite-volume partition function one gets (z/^) = VEm/Nf and hence the 
linearly raising behaviour (^ of the topological susceptibility, valid for x^l. Analogously, 
one may compute x from in the other two regimes, since Leutwyler and Smilga give, in 
their landmark paper |^, also the generalization of ([TBI) to arbitrary x: 



7XPT 



Iu{x) 

h-lix) 



Iu+l{x) 

Ux) 



Iu-Nf+l{x) Ii,^Nf+2{x) 



Iu+Nf~l{x) 
Iu+Nf~2{x) 

Iu{x) 



(21) 



It is easy to check that, upon using the asymptotic form /,^(x)^e^(l — (z/^— l)/(8x))/-\/2vrx 
for X ^ 1 of the modified Bessel function of integer order, the representation (21) of the 
partition function reduces to (|16[) for x^l, hence reproducing the linearly raising form (^ 
of X ill the large LS-regime. Furthermore, from the representation Iu{x) ~a;''''/(2''^'r(|i/|+l)) 
appropriate in the opposite case x<^ \ one gets the asymptotic behaviour [H 



N, 

X cc m f 



(x<l) 



(22) 



characteristic for the small LS-regime. For an x value which is neither small nor large 
compared to one, a numerical evaluation of (z/^) = Sz/^Z^^"'"/SZ^^"'" is needed, and the 
result may be compared to the asymptotic behaviour (|^) or ( p2t) for x>l or a;< 1. One of 
these comparisons is demonstrated on the l.h.s. of Fig. 6. There are two important lessons 
we can learn from that graph. First, finite-volume effects have a decreasing effect on the 
topological susceptibility. Second, the topological susceptibility as derived from ( pl| ) and 
the large-x asymptotic form (^ are reasonably close for moderate x-values already. To be 
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In addition, bosonic/fermionic fields are supposed to be periodic/antiperiodic in that direction. 
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Figure 6: LHS: Topological susceptibility (without instanton effects) [in MeV^j versus quark mass 
[in MeVj as derived from (full curve), and the large-x asymptotic form (dashed line) in 
a box of V = 2(1.5fm)'* with T, = (288 MeV)^. Since Nf = 2, the small-x asymptotic form is 
a parabola (cf. RHS: Illustration of the scenario in which the topological susceptibility 

in a finite volume (full) "overshoots" the infinite-volume curve (dashed) for some intermediate 
quark masses, hence explaining the CP-PACS "flatness problem" showing up in Fig. 3. In this 
illustration, the "granularity" of the QCD vacuum, as discussed in the text, is already taken into 
account. This "overshooting" scenario is excluded, since Enite volume effects may only have a 



decreasing effect on x (cf- LHS and the "granular" construction (27) discussed in the text) 



definite: For V = 10. 125 W and S = (288 MeV)^, the mark X = 1 is reached at a quark 
mass m = 6.26MeV, and there the finite-a; value is about 50% smaller than the asymptotic 
prediction. Taking the quark mass 10 times larger, m = 62.6MeV (i.e. x = 10), the finite-a; 
value lies less than 10% below the asymptotic (large-x) prediction (cf. Fig. 6). 

The first statement is important, because it rules out a possible explanation of the 
"fiatness problem" in the CP-PACS determination of the topological susceptibility on the 
smaller (16'^x32, P = 1.95) lattice, apparent in Fig. 3: Knowing that the topological suscepti- 
bility drops, in the small LS-regime, with a power-like behaviour - cf. (^2]) - and thus stays 
massively below the linear curve (^), one might have suspected that x could "overshoot" 
and lie above the linear curve in an intermediate regime. Such a hypothetical behaviour, as 
illustrated in the r.h.s. of Fig. 6, which would provide a simple explanation why the leftmost 
point in Fig. 3 is so high, is definitely ruled out by the graph on the l.h.s. of Fig. 6; in other 
words, the "chiral constraint" indicated in Figs. 3 and 4 is to be taken seriously. 

The second statement is important, because it implies that in a lattice study the mini- 
mum box length needed to avoid finite-volume effects in the topological susceptibility is 
much smaller than what one might have suspected by analogy with the standard pion cor- 
relation length criterion. The standard criterion says that in the broken phase the minimum 
box length L^in has to be larger than by up to a factor 4 to avoid finite-volume effects 
in hadronic Green's functions. If the same factor (per linear dimension) would be relevant 
for the LS-parameter x, the latter would have to be of order 4^ = 256, which, as one can 
see from Tables 1-5, is hardly the case in any of the simulations discussed in this article. 
Hence from the l.h.s. of Fig. 6 one takes the good news that the actual value needed is much 
smaller; already x > 10 seems sufficient to make the curve lie reasonably close to the large-x 
line. Furthermore, comparing x as predicted by the QCD partition functions (|16D and (p]) 
offers the potential benefit to quantify finite-volume effects in the topological susceptibility. 

The next step is to recall that there is an important limitation to either topological 
susceptibility curve shown in the l.h.s. of Fig. 6, i.e. this limitation is inherent to both 
the solid curve (numerically obtained from ([2T|)) and the dashed line (derived from (p!6|)): 
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Either function grows without bound, i.e. they both violate the other asymptotic constraint, 
X ~^ Xoo for m — i> oo, as discussed in Sec. 2 and illustrated in Fig. 1. 

In such a situation it is crucial to reflect on the physical basis of the second constraint, 
valid at large quark masses. It is easy to see that the relevant physics is due to instantons; 
they are sufficient to get an understanding of the numerical value of the quenched topological 
susceptibility, Xqo- Assume the total four- volume V of the box to be sufficiently large, so 
that it can be divided into N hypercubic boxes of volume Vq each, and - without loss of 
generality - even. In the Instanton Liquid Model (ILM) the typical volume occupied by 
an instanton or anti- instanton is — 1 fm ||T^, so we expect Vq to be roughly of that size. 
Assume now that each one of these boxes is populated, at random, by either an instanton 
or an anti-instanton. What is the resulting topological susceptibility ? The answer is easy: 
The topological charge distribution is a modified binomial distribution (note that z//2 G Z) 



= (23) 




which, after it has been smoothed to fill the odd-u entries too, may be approximated by 

^ exp{-^} , a' = N = ^ (iV»l). (24) 



2(t2^ ' Vo 
This distribution, in turn, implies (z/^)~V/Vo, and henceforth 

Xoo ^ ^ , (25) 

which means that, if the cell volume Vq was identical to V^nst = 1 fm"^, the quenched topolo- 
gical susceptibility would evaluate to 1.52 lO^MeV^. Hence, all we need to do, in order to 
get a more standard value of the quenched topological susceptibility, is to increase the size 
of the little boxes, i.e. to pack the instantons a little bit looser: 

Vo ~ 1.44 fm^ ~ (1.10 fm)^ Xoo ^ 1-05 lO^MeV^ ~ (180 MeV)^ . (26) 

The basic lesson from this exercise is that it is possible to get a rough quantitative under- 
standing of the quenched topological susceptibihty from a simple statistics consideration 
involving nothing but the parameter Vq, the average volume occupied by an instanton, hence 



omitting all the dynamics between instantons which is incorporated in the ILM [|18 



It is now easy see why the topological susceptibility as derived from the finite-volume 
partition function (^) is, in general, too large. To this end it is useful to consider the 
distribution of topological charges as predicted by ( PT| ) in a box volume typical for full 
QCD simulations, say \^ = 2 (1.5 fm)"^ = 10.125 fm"^, at a quark mass of, say, 250 MeV. As 
one can see from Fig. 7, in this regime the partition function Z^^"^ as given in ( pi]) is 
wider than the partition function Z^^^ as given in (p^. In particular, for sufficiently high 
\v\ the former indicates a sizable population of that topological sector, whereas the latter 
is almost zero. This means that the (leading order) chiral prediction neglects the basic 
fact that instantons form an incompressible fluid ||18|], i.e. it is impossible to pack such 
a large number of instantons into the box as the tail of the Z^^"^ distribution suggests. 
This calls for an improvement on Z^^"^ which incorporates the unique feature of Z^^^ that 
instantons occupy a fixed volume. In the absence of an exact analytic solution it is natural 



to assume that the two suppression mechanisms (captured in (pT| , p^), respectively) operate 
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Figure 7: Normalized QCD Bnite-volume partition function Zy as predicted by the Leutwyler- 
Smilga form ( |2j|j (top left) and by the instanton inspired form ( [2^ ) (top right) for V = 2 (1.5 fm)^ 
and m = 250MeV. Following the XPT prediction, the topological sectors with >8 should be 
populated, but the ILM statistics consideration demands that they are heavily suppressed, since 
instantons cannot be compressed. The construction (27) assumes that the two high-li'l suppression 
mechanisms do not interfere. The resulting "QCD finite-volume partition function with granula- 
rity" is narrower than either one of its predecessors (bottom left), and the topological susceptibility 
curve derived from it (bottom right) shows all essential features: power-like behaviour for x <1, 
almost linear behaviour for x>l and I^m/Nf <Xoo, saturation for T,m/Nf>Xoo- 



independent of each other, i.e. to construct the QCD finite-volume partition function by 
just multiplying^ the two asymptotic^ partition functions 

^QCD ^ ^XPT . ^ILM ^ (27) 

as is illustrated in the lower left plot in Fig. 7. As one can see, the resulting distribution is 
narrower than either one of the asymptotic distributions it was constructed from. This still 
holds true in the opposite regime of small quark masses where the relationship is reversed 

^^It is worth mentioning that this is reasonable regardless of what are the respective energy and entropy 
components in the two partition functions. For a single degree of freedom (|27| ) means that the energy 
differences w.r.t. Eq add up. For a larger number of degrees of freedom either partition function is, in 
general, the product of a Boltzmann-typc energy suppression and an entropy factor. In our case Z^^^ 
is likely of this "mixed" type (even though the axial WT-identity from which it takes its origin does not 



disclose any information on this point), but the second factor Z^^ as constructed in (24) is pure entropy. 

'^'^Note that the multiplicative rule ( p7| ) w.r.t. v is equivalent to a folding prescription in the basis w.r.t. the 
dual variable 9: The likelyhood that the system is found in the interval [9^^^, 6'inax] is Z^^'^{e') ■ Z^^^{e") 
integrated over all 9' , 9" with the constraint that the sum 9' + 9" is in the interval [6*111111, 6'max]- 

^^The meaning is: appropriate, by themselves, for to — > or m — > 00, respectively. The choice of the 
representation (^ ) for Z^^"^ and of (|2j) for Zj}"^ makes sure that (|7|) is appropriate for arbitrary x and 
any V large enough to accommodate several instantons (e.g. V> lOfni^). 
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(there Z'^^'^ is narrower than Z^^^ which is quark mass independent). As a result of this 
the topological susceptibility is smaller than both the XPT-prediction and the ILM-inspired 
value. This is particularly easy to see in the limit of large box volumes; there either partition 
function entering the r.h.s. of (|27|) tends to be gaussian (with variance a"^ = VT,m/Nf and 
<^^ = ^Xoo, respectively), and the multiplication prescription (pT]) then implies 



xiNf, m) = l/(iV//Sm + 1/xoo) , (28) 

which indeed fulfills xi^f-, ^) ^ min(Sm/A'/, Xoo)- Note that ( pSj) is nothing but the member 
n=l out of the phenomenological family (^. In addition, two points should be mentioned: 
First: The nice feature that finite-volume effects in x disappear for moderate (x, V") combi- 
nations already (say for x > 10, > 10 fm^) has survived. This follows from the construction 
(P7|) of the "granular" QCD finite-volume partition function, but one may also verify it by 
carefully comparing the last graph in Fig. 7 to the dashed ( "infinite volume" ) curve on the 
r.h.s. of Fig. 6 - the graph in Fig. 7 was generated for V = 2 (1.5 fm)'* = 10.125 fm^, hence 
reaching x = 10 at m~62.6MeV. Second: As a specific finite- volume effect the topological 
susceptibility curve may show, in a certain (typically narrow) window of (typically small) 
quark masses, a positive curvature - a feature absent in the infinite- volume limit. 

We shall conclude this paragraph with a short reflection on what we have achieved with 
the "granular" construction ( P?] ) of the (true) QCD finite-volume partition function. 

On the conceptual level a specific observable, the topological susceptibility x, has been 
linked to a more general concept, the QCD finite- volume partition function with granularity. 
The name reminds one that the expression (^) which resulted from a careful evaluation, at 
finite x, of the QCD partition function in terms of the chiral Lagrangian to order 0{p^) p| 
got modified to account for the "granular structure" of the QCD vacuum, due to instantons. 
Note that the "granular" partition function depends separately on V and m - not only on 
the product of them, as (|T^) and (^) did. The construction is not ad hoc, because the 
distribution Z^^^ predicts, besides another observable, namely the mass dependence of 
the sectoral chiral condensate in a finite volume, hence promising an interesting generaliza- 



tion of the results by Damgaard |T9[. On the other hand, it is clear that the construction 
is neither a model nor exact, since subleading effects (where the dynamics of QCD enters) 
have been neglected; the instanton effect taken into account is pure entropy^ 

On the practical level, the construction (^) has managed to select, for purely theoretical 
reasons, out of the phenomenological family (^) the specific member (^) as the one which 
is, in the limit of large box volumes, most promising. This provides us with an overall fitting 
curve with just two parameters which is supposed to be reasonable for arbitrary quark mass 
values. Furthermore, at any given box volume a detailed comparison of the topological 
susceptibility curve as predicted by ( P?] ) to the "infinite volume" form (^) allows one to 
come up with a quantitative assessment of finite-volume effects on x- This may prove useful 
in future simulations of full QCD, since it enables one to choose the box volume just as 
large as needed in order to keep finite- volume effects tolerable, hence providing some help 
in the challenging task to allocate computational resources as reasonably as possible. The 
limiting factor is, of course, that even the demonstrated smallness (within the accuracy of 
the construction (|27|) ) of finite- volume effects in x is no proof that finite- volume effects in 
other interesting observables would be small. 



^^In this respect it should be noted that the suffix "ILM" is, strictly speaking, a misuse of the name, 
since the real instanton liquid model does account for the dynamics between instantons [nSl . 
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6 Fits, Cuts and Speculations 



We are now ready for a detailed analysis of the available data for the topological suscepti- 
bility in 2-flavour QCD. 

The most valuable output of the theoretical considerations presented above is a well- 
motivated functional form against which we may fit the data over the entire range of quark 
masses. This is particularly useful because, as mentioned in the second section, the bulk 
of the data is in the transition regime; hence functional forms which are inspired by an 
expansion appropriate in one asymptotic regime are not guaranteed to yield reliable results 
for the physical parameters we are interested in. 

For large enough box volumes (for details see below) formula ( pH]) will be applicable. It 
contains two parameters, S and Xoo, on which the simulations will disclose some information. 
Regarding the chiral condensate in the chiral limit, it is important to note that the way 
it will be determined (via fitting (^Sf ) to the lattice data) is logically independent from the 
usual approach of measuring the condensate at various quark masses via the trace of the 
Green's function of the Dirac operator and then extrapolating to zero mass. The procedure 
to determine S from the topological susceptibility curve does not require a subtraction or 
difficult-to-know renormalization constants that arise in more conventional approaches. 

Fig. 8 displays the data by CP-PACS (on the 24^x48 lattice only), UKQCD, SESAM 
and the Pisa group together with individual interpolating curves which represent both 
reasonable fits of the data to the functional form (28) and of the inverse data versus the 
inverse argument to a straight line. The respective values for E and Xoo are tabulated in 
Table 6 (along with further details on how they are determined). A naive average (performed 
in the top line in either part of the table) yields E ~ (425 MeV)^ and Xoo — (203 MeV)'^, 
respectively, or E ~ (395 MeV)'^ and Xoo — (210 MeV)^ if the values suggested by the Pisa 
group study are not included. A direct fit to the combined data (see Fig. 9) suggests 
E ~ (484 MeV)^ and Xoo — (194 MeV)^. It seems that the outcome of this combined analysis 
is influenced, to a huge extent, by the Pisa data - even stronger than the naive average. 

To ensure the quality of the final phenomenological parameters it seems necessary to 
require the data points to meet certain criteria to be included in the analysis. What we have 
learned from the previous two sections is that we should not only require a to be reasonably 
small to keep lattice artefacts under control, but we should also demand the LS-parameter 
X and the volume V to be reasonably large to keep finite-size effects on x moderate and 
to be able to accommodate a reasonable number of instantons. What I dare proposing, for 





CP-PACS 


UKQCD 


SESAM 


PISA 


3-avg 


4-avg 


all 


sel 




24.6 


46.4 


16.4 


5.6 


29.1 


23.3 


15.8 


24.1 


[MeV] 


162. 


118. 


199. 


340. 


149. 


167. 


203. 


164. 


EV3 [MeV] 


418. 


338. 


478. 


684. 


395. 


425. 


484. 


420. 




0.063 


0.083 


0.076 


0.035 


0.074 


0.064 


0.054 


0.063 


Xoo [109 Me V^] 


1.7 


2.2 


2.0 


0.9 


1.9 


1.7 


1.4 


1.7 


X]L' [MeV] 


202. 


216. 


211. 


174. 


210. 


203. 


194. 


202. 



Table 6: Results of the fits of the data by the individual groups (Fig. 8) or of the combination 
of all or a selection of them (Fig. 9) to the functional form (p^j and of l/(xro) as a function 
of l/{M^rQ) to a straight line, weighting the results 2:1. For both the inverse chiral symmetry 
breaking parameter (upper part) and the quenched topological susceptibility (lower part) the 
average over the first 3 or 4 collaborative results is taken in the first line. Successive lines are 
derived using ro = 0.49fm and, to convert F.,^ into S, the values for M-,j,m quoted in footnote 5. 
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Figure 9: Fitting all data (top) or only those which pass the cuts ([2^) (bottom) to the functional 



form (28). Either way S seems to he much larger than expected, as one can see from comparing 
the fitted curve to the one anticipated - based on S = (288 MeV)^ and the most optimistic value 
Xoo = (200 MeV)^ - which is included in the second plot. Data taken from ^, ||, |^ |^ |^. 



these reasons, is to apply the following cuts on the data 

(a) a<0.15fm , (b) x>10 , (c) V>10im\ (29) 

irrespective of the circumstances^. From Tables 1-5 one sees that this battery of cuts 
eliminates the CP-PACS data on the smaller grid and those by the Pisa group. In addition, 
the UKQCD data set looses its most chiral point, and for the SESAM data only those 
generated on the larger grid and the two least chiral ones on the smaller grid survive. Table 
6 contains, in its last column, the results of a direct analysis of the combined ("selected") 
data which have survived the cuts (^). The net outcome almost coincides with the naive 
average of the results by the individual groups; in the light of the sizable deviations^ 

^^Cut (a) is unfair w.r.t. the CP-PACS data; it does not care whether the gluons are improved or not. 
^^This is the reason why no error bars are given in Tables 1-6. It seems that overall systematic effects 
dominate over statistical fluctuations. 
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between the results by the various groups (cf. Table 6) I quote 

S = (390...430MeV)^ or = 145...170 MeV (30) 

Xoo = (200...210MeV)^ (31) 

as the current result for S and Xoo when determined via the topological susceptibility curve 
(and without continuum extrapolation) in 2-flavour QCD. 

For the quenched topological susceptibility the result (^) is slightly larger than Xoo as 
determined in previous direct studies in the SU{3) theory [^. This is not disturbing, since 
the generally accepted value Xoo — (190...200 MeV)^ is after a continuum extrapolation, 
whereas the value (0) is gained at finite a; hence the difference may - as it is "small" - be 
attributed to the 0{a^) term showing up on the r.h.s. of eqn. (|r^). 

For the chiral symmetry breaking parameter the situation seems less pleasing; the result 
(^) is substantially larger than the generally accepted ~ 93MeV, S ~ (288 MeV)^, as 
determined experimentally or via the GOR- relation, respectively (cf. footnote 5). This 
excess may either be due to lattice artefacts or it may reflect an interesting A^-dependence 
of QCD close to the chiral limit. We shall discuss these two possibilities in due course. 



Regarding lattice artefacts as a possible reason why (|30D exceeds our expectations, little 
can be said beyond the points already raised in section 4. The issue is not whether lattice 
artefacts are present - they are, since no continuum extrapolation has been attempted - 
but rather whether they are numerically under control, i.e. whether the final result (|30D is 
largely influenced by discretization effects or not. The point is that we have just seen that 
the value ( pT]) for Xoo happened to be in rough agreement with our expectations even without 
a continuum extrapolation, and if lattice artefacts are "small" for the quenched topological 
susceptibility, it is natural to expect their influence on the other measured quantity, S (or 
-Ftt) to be tolerable too, and this is apparently not the case. 

Lattice artefacts have been categorized, in section 4, as (a) scaling violations and (b) 
chirality violation effects. While these two types of lattice artefacts may, in general, be inde- 
pendent of each other, with Wilson fermions and its descendents (as used in the simulations 
discussed here) they go together, i.e. they jointly disappear in the continuum-limit and only 
then. Since in all the data which have passed the cuts (^) smaller quark masses go, in 
tendency, together with smaller lattice spacings (except for the UKQCD data which are par- 
tially matched, cf. Table 3), we would expect the deviations from the dotted curve (which 
reflects our expectations based on (^) together with S = (288 MeV)^, Xoo = (200 MeV)"^) to 
become smaller for low quark masses, if (a) or (b) (or both of them) represent the source of 
the supposed excess. The bottom part of Fig. 9 (where finite-volume effects are negligible) 
may contain a hint in this direction: If we suppress the points by the SESAM collaboration 
(which have huge error bars anyway) the deviation, by absolute standards, of the data from 
the dotted curve seems to increase towards the r.h.s. (where a is, in tendency, larger). In 
addition, it looks like there is a slight tendency for the UKQCD data to lie beneath the 
CP-PACS-2 data. This nicely corresponds with the order of typical lattice spacings as used 
in these simulations (cf. Tables 2-5), even though unequal efforts have been spent on im- 
proving the action. These observations certainly support the view of lattice artefacts as the 
primary reason why the value (150) for the chiral order parameter exceeds our expectations 
by roughly a factor 2 (which means a factor ^/2 for F^,), but they are far from conclusive. 



The second possibihty is that the values (|30| , pl| ) are indeed correct, in other words that 



the problem is not with the lattice determinations of S (or F^r) and Xoo, but rather with our 
expectations - specifically with our knowledge regarding the chiral order parameter. This 
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might come as a surprise, since it is generally believed hat the value F^^ — 93 MeV or S ~ 
(288 MeV)^ is rather well known experimentally or through the GOR-relation, respectively. 
The point is that the quoted (largely different) values need not necessarily be in conffict 
with each other, since (|30|) is determined in Nf = 2 QCD, while the experimental value is 
in QCD with0 Nf = 3. The statement is that there might be an iV/- dependence in the 
topological susceptibility near the chiral limit beyond the one indicated in (^). In other 
words, there might be an A^j- dependence in the chiral condensate itself which makes S2 as 
determined in the lattice studies under discussion differ from the experimentally "observed" 
S3, contrary to the standard belief that SjVj is0 largely independent of A^. 

The idea that S might depend on Nf is neither new nor exotic. It is known that QCD 
has a "conformal window" with restored chiral symmetry if the number Nf of massless or 
light flavours is appropriately chosen, i.e. Nf^^^^<Nf<Nf^^^^ [^. While the upper critical 



number for this phase is the number of flavours where QCD looses asymptotic freedom, i.e. 
^/crit = 33/2 (for A^ = 3) [^, the lower critical number is much under debate. Originally, 
the Orsay group argued that iV^'crit could be so low (say A^'crit — ^(3), as opposed to the 
more standard view Aj^rit — 0{7) [^]) that the chiral condensate in ordinary 3-flavour 
QCD might be marginal or even vanishing, hence necessitating the extension of the chiral 
Lagrangian approach |]IU| to the frameworkp^ of "generalized chiral perturbation theory" 



(GXPT) Meanwhile, they accept that the lower critical bound might be somewhat 

larger, say A^crit = ^(5), but the interesting question remains whether there is any trace 
of the adjacent "conformal window" for lower A(/- values where QCD is in the conflned 
phase and the chiral symmetry is broken - the latter most likely signaled by the lowest 
dimensional possible order parameter already, i.e. though S >0 (distinctively). If in this 
"ordinary" phase S would indeed get substantially reduced when Nf is increased by one 
unit, such a behaviour could be interpreted as a hint that the real world with Aj = 3 is 
already "close" to the lower bound A^'crit; as argued by the promoters of GXPT p4|. Hence, 
if the values quoted in ( PU] ) would indeed represent the ultimate truth for Nf = 2 (while still 
~ 93 MeV or S ~ (288 MeV)^ for Nf = 3), this could be taken as an indication that 
the real world is not too far away from the lower end of the "conformal window". The 
important point is, of course, that such a hypothesis is accessible to direct tests on the 
lattice^, and due to those the current body of evidence in favour of a relatively low A^'crit 
p6| , ^ will either be whiped out or corroborated^. As an ironic twist one should note that 
the value (|30|) which we supposed as "ultimate" in this scenario and which pushed us into 
these speculative thoughts is, by its numerical size, anything but marginal. 



^^This is not to be confused with the fact that the (non-strange) M^r would hardly change if the strange 
quark would be made infinitely heavy. The statement is that -FV might feel the infiuence of the strange 
sea-quarks, hence spoiling the generally believed "stand alone" quality of the sector with 2 light flavours. 

^^The understanding is: for the range of flavours we are considering here, i.e. for Nf = 2..A. 

^°With Nf close to Nj^^^^ the supposed order parameter S becomes small, and the "correction terms" on 
the r.h.s. of the GOR-relation M^F^ = 2I]m -I- ... turn larger than the "leading" term. As a result of this, 
the chiral expansion must be reordered - at the price of reducing the predictivity of the theory at a given 
order in the expansion in external momenta. For a review of GXPT see . 

^^A simple test would be to measure xl"^) in QCD with Nf = A degenerate flavours and to show that the 
topological susceptibility curve has, for m 0, an asymptotic slope which is less than half of that in the 
bottom part of Fig. 9. The ultimate test is, of course, to show that ~ despite the supposed correctness of 
® " ^^^^ fV~93MeV, I]~(288MeV)3 in 2-hl-flavour QCD. 

^^^Intuitively, one might think that the result ( pO| ) of the present analysis is strong evidence against this 
scenario, but things are not that simple: If A^'^rit turns out so low that even E2 is zero or marginal (i.e. 
not sufficiently large that the first term on the r.h.s. of the GOR-relation dominates; cf. footnote 20), then 
the result (pO[) itself is affected, since in such a case E may not simply be read off from x(m). 
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7 Summary and Outlook 



In this paper I have confronted recent lattice data for the topological susceptibility x i^i 
Nf = 2 QCD as a function of the (dynamical) quark mass with theoretical expectations. 

The ultimate goal of these studies is to understand how QCD manages to reduce the 
prominent ^-dependence of quenched QCD up to the point where it is eliminated, if the 
chiral limit is taken. Dynamical simulations, by themselves, provide little insight into the 
mechanism behind this fascinating behaviour; what they can do, however, is to document 
that the topological susceptibility gets suppressed (and hence the 6'-dependence reduced) 
if the dynamical quark masses are tuned sufficiently small. Pure theoretical thought, on 
the other hand, may provide interesting bounds and it may yield information about some 
asymptotic behaviour, but typically the information produced is qualitative, i.e. the relevant 
constants have to be "borrowed" from elsewhere. 

The approach taken here is to combine the lattice data with analytical insight. This is 
complicated by the fact that the preferred ranges of quark masses are quite different. To 
be on safe analytical grounds one likes the quarks to be either sufficiently light so that the 
flavour- singlet axial WT-identity enforces the linear relationship (^ or sufficiently heavy so 
that the topological susceptibility effectively coincides with its quenched counterpart. The 
lattice data, on the other hand, happen to be precisely in the intermediate regime which 
- based on our preliminary knowledge on S and Xoo in QCD - was estimated to be at quark 
masses of order 88 MeV in the 2-fiavour theory and where the topological susceptibility 
follows neither one of the two asymptotic patterns. 

The key idea on which this article is based is that it is possible to combine the leading 
physics resources of either regime without invoking specific model assumptions. This is 
done on the more general level of the finite-volume QCD partition function Z„ from which 
the topological susceptibility may be derived. What enters the construction ( P7| ) is the 
assumption that the two mechanisms which suppress configurations with higher w.r.t 
those with lower |z/| in the small- and the large-mass regime, respectively, operate indepen- 
dently of each other. For small enough quark masses the relevant mechanism is captured 
in the leading order chiral Lagrangian (which in turn implements the constraint imposed 
by the axial WT-identity), and the suppression gets manifest upon evaluating the func- 
tional integral over the pseudo-Goldstone boson manifold. For large enough quark masses 
instantons are the relevant degrees of freedom, and an argument is pushed forward that 
the suppression of higher topological sectors follows from a pure entropy consideration, if 
instantons are taken incompressible. To the extent that instanton interactions may be ne- 
glected and the relevant physics is, on that side, indeed pure entropy, the hypothesis that 
the two suppression mechanisms are independent is then correct by construction. 

On the practical level, the outcome of this analysis is twofold. First, for large enough box 
volumes the functional form (^) is established as a well-motivated fitting curve, suitable 
for the entire range of quark masses encountered in present and future studies of QCD 
with Nf>2 (and Nf < Nj^^^^). It contains only two fitting parameters, and both of them 
have a direct interpretation in terms of (semi-)phenomenological quantities (S and Xoo)- 
Second, the two different choices for the chiral finite-volume partition function Z^^'^ (i.e. 
(0) vs. (0)) entering the construction of Z^'^^ allow one to quantitatively assess - within 
the validity of (|27|) - the impact of finite-volume effects on the topological susceptibility. 
The somewhat surprising outcome is that they are smaller than present days statistical 
uncertainties for relatively small values of the LS-parameter x already, say for a;>10. 

Making use of the functional form (|2|) the data by CP-PACS, UKQCD, SESAM/TXL 
and the Pisa group have been analyzed. The results for the asymptotic parameters S and Xoo 
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seem to vary more than what one would have expected from the error bars of the individual 
data points. While the quenched topological susceptibility Xoo is in loose agreement with 
both the estimate from the Witten-Veneziano relation and previous direct determinations 
in the quenched theory, the result for the low-energy constant S = — limm^o limv^oo('^'?/^) 
outstrips the expected value by a factor 2. Here the expectation relies on the GOR-relation 
and the lattice determination of the non-strange quark mass; would one use a more tradi- 
tional estimate of (m„ + mrf)/2, the factor by which S as determined from the topological 
susceptibility curve exceeds our expectations would be larger. Even restricting the set to 
contain nothing but "high quality" data (which have passed certain cuts to ensure that 
effects due to finite lattice spacing and finite box volume should be small) hardly changes 
the result of the analysis - the Nf-io\d asymptotic slope for Nf = 2, E2, stays high. 

Two possible reasons for this excess have been discussed. The simple explanation (sup- 
ported by some inherent features of the data) is that it is due to lattice artefacts. The 
alternative one (which is mainly discussed because of it conceptual appeal and its potential 
to initiate dedicated studies) is that S2 is indeed much larger than the GOR-relation sug- 
gests, indicating a larger-than-believed A^-dependence of low-energy QCD. In either case, 
it seems fair to say that at least in 2-fiavour QCD chiral symmetry is predominantly broken 
through a (distinctively nonzero) chiral condensate 



The implication is, of course, that further studies are needed to clarify the situation (cf. 
From the analysis one gathers that there are basically three reasonable strategies for 
going more chiral: to simulate (z) at fixed (3, [ii) at fixed a, {in) at fixed x, whereas the 
possibility {iv) to simulate at fixed rh seems to yield more questionable results, since in this 
approach (3 gets so much relaxed that the lattice gets intolerably coarse for the most chiral 
points. While it is clear that in principle (i.e. if a continuum extrapolation were possible) 
any of these possibilities would yield the same physical answer, in practice the choice may 
make a sizable difference, and it is important to keep in mind that the "optimum" choice 
may not only depend on the regime of quark masses and box volumes one wishes to cover, 
but also on the actions used. In this respect it is clear that - even if one cannot afford a 
continuum extrapolation for all the data points constituting the x = xij^) curve - a scaling 
study at least for a single quark mass value might help to get a more quantitative assessment 
of discretization effects and also to indicate which might be the most promising combination 
of actions in the gluon and fermion sectors, respectively. If progress turns out to be slow 
in the Nj = 2 system, studying QCD with Nf = 4 might be an interesting digression - not 
only because S4 would tell us whether the "conformal window" is close, but mainly for 
the practical reason that the "transition regime" occurs at higher quark mass values and is 
therefore cheaper to study. What we are ultimately interested in, however, are the values 
gained in a more ambitious project, namely the simulation of 2-|-l-flavour QCD. 
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